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ABSTRACT 

^ The population of rich galaxy clusters evolves much more rapidly in a universe 

with critical density than in a universe with low density. Thus, counts of clusters 
\^ ' at intermediate redshift offer the possibility of determining the cosniological density 

parameter, fio, with a minimum of assumptions. We quantify this evolution using 
the Press-Schechter formalism which we extend to flat cosmological models with a 
cosmological constant, Ag = 1 — ilo- Using new large N-body simulations, we verify 
that this formalism accurately predicts the abundance of rich clusters as a function 
of redshift in various cosmologies. We normalise the models by comparing them to 
the local abundance of clusters as a function of their X-ray temperature which we 
rederive from data compiled by Henry & Arnaud. The resulting values of the rms 
O ; density fluctuation in spheres of radius Sh-^Mpc are cts = (0.50 ± 0.04)f)(7°"'^+°-^°"" 

^ ! if Ao = and erg = (0.50 ± 0.04)f7o if Aq = 1 - VLq. These values depend 

only weakly, and almost not at all if f2o = 1, on the shape of the power spectrum. 
fH , We then examine how the distributions of mass, X-ray temperature and Sunyaev- 

Q^' Zel'dovich decrement evolve as a function of We present the expected distributions 

at z = 0.33 and z = 0.5 and the predicted number counts of the largest clusters, both 
in space and in projection on the sky. We find that even at z — 0.33, these distributions 
depend very strongly on JIq and only weakly on Aq. For example, at this redshift, we 
^ \ expect 20 times as many clusters per comoving volume with M > 3.5 x lO^^h"^M0 

and 5 times as many clusters with kT > 5 keV if JIq = 0.3 than if fio = 1. The 
splitting in the integrated counts is enhanced by the larger volume element in low fio 
, models. There is therefore a real prospect of estimating f^o from forthcoming surveys 

j_j ■ of intermediate redshift clusters that will determine their masses, X-ray temperatures 

' or Sunyaev-Zel'dovich decrements. 
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1 INTRODUCTION 

Galaxy clusters are exceptionally useful tools for estimating 
fundamental cosmological parameters. Their utility stems 
largely from their relative dynamical youth. In hierarchical 
clustering theories, clusters form by the gravitational ampli- 
fication of primordial density fluctuations, usually assumed 
to have an initially Gaussian distribution of amplitudes. 
Rich clusters correspond to the rarest collapsed objects in 
this distribution and their abundance varies very rapidly 
with properties such as mass or potential well depth. 

The mass within the virial radius of a rich cluster 
is very close to the mass enclosed within 
a sphere of radius 8h~^Mpc in the unperturbed universe 
(Evrard 1989, White, Efstathiou & Frenk 1993). Because 



of this, the present day abundance of rich clusters directly 
reflects the amplitude of density fluctuations on a scale of 
~ 8h~^Mpc and can be used to measure this amplitude 
with a minimum of assumptions. In general, this measure 
depends only weakly on the amplitude of fluctuations on 
other scales. It does, however, depend on the value of the 
cosmological density parameter, fio- Thus, the observed lo- 
cal cluster abundance flxes the value of (Tg, the rms density 
fluctuation in spheres of radius Sh^^Mpc, as a function of 
no. 

The temporal evolution of the cluster abundance is de- 
termined by the rate at which density perturbations grow. 
This, in turn, depends primarily on the value of Q, and, to 
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a lesser extent, the value of the cosmological constant, 
and the shape of the power spectrum of density fluctuations. 
In a low density universe, fluctuations cease to grow after a 
redshift z ~ {^o^ — (e.g. Peebles 1980), resulting in a 
cluster population that evolves very slowly at low redshift. 
In an fio = 1 universe, on the other hand, density fluctu- 
ations continue to grow even at the present epoch and so 
the cluster population is still evolving rapidly. Measurement 
of the cluster abundance at moderate redshifts provides a 
powerful method for estimating SIq. 

In order to apply these tests, it is necessary to pre- 
dict and measure the abundance of clusters as a function of 
some property such as mass. X-ray luminosity or the equi- 
librium temperature of the X-ray emitting intracluster gas. 
Theoretical predictions and observational measurements are 
subject to different uncertainties. The cluster mass within 
the virial radius is the simplest quantity to predict theoreti- 
cally, but one of the hardest to measure reliably. Traditional 
virial analyses are prone to contamination by projection ef- 
fects (Frenk et al. 1990; Van Haarlem, Frenk & White 1996; 
Mazure et al. 1995) and X-ray data do not generally extend 
to the virial radius. A novel and highly promising technique 
to measure cluster masses employs the shape distortions of 
background galaxies produced by weak gravitational lens- 
ing in the cluster potential (Kaiser & Squires 1993, Fahlman 
et al. 1994, Wilson, Cole & Frenk 1996, Seitz & Schneider 
1995 and references therein.) The limited field of view of the 
current generation of CCD cameras, however, restricts such 
measurements to the inner parts of clusters. 

The simplest quantity to measure empirically is the 
cluster X-ray luminosity. However, since the bremsstrahlung 
emissivity per unit volume is proportional to the square 
of the gas density, the total power radiated is very sensi- 
tively dependent upon the distribution of gas in the cluster 
core. This is difficult to model, particularly at high redshift 
(Evrard 1990; Navarro, Frenk & White 1995 and references 
therein). Nevertheless, X-ray luminosity provides a conve- 
nient means to select complete samples of galaxy clusters. 
In contrast, the temperature of the intracluster gas can be 
predicted quite reliably using modern hydrodynamic simu- 
lations. These show that as a cluster collapses, the gas is 
shock heated to the virial temperature and rapidly settles 
into hydrostatic equilibrium with an approximate isothermal 
structure (Evrard 1990, Cen & Ostriker 1994, Bryan et al. 
1994, Navarro et al. 1995). Average (luminosity- weighted) 
X-ray temperatures have now been measured for fairl y large 
sample s of clust ers, with the Eins tein Observ atory (Henry 
& Arn4ud 1991]), and EXOSAT ([Edge 1988|), and the ra- 



dial variatio n of temperature is now beginning to be probed 
with ASCA ( Markevitch et al. 1996 ). Early results show that 
the gas is indeed approximately isothermal. Finally, another 
observable that is also insensitive to the detailed distribu- 
tion of the gas within the cluster is the AT/T decrement in 
the cosmic microwave background radiation produced by the 
Sunyaev-Zel'dovich (S-Z) effect (Sunyaev & Zeldovich 1972). 
The line-of-sight decrement from the cluster depends only on 
the product of the column density of gas and the temper- 
ature and is independent of the cluster redshift. Sunyaev- 



* We express A in units of 3Hg, where Hq = lOOhkms ^Mpc ^ 
is the present value of the Hubble constant 



Zel'dovich decrements have now been reliably measured for 
several clusters (eg Grainge et al. (1993), Wilbanks et al. 

(1994) and Birkinshaw & Hughes (1994)). 

In this paper, we exploit those cluster properties which 
are best determined theoretically and observationally in 
order to constrain the values of cosmological parameters. 
We first use the distribution of X-ray temperatures at the 
present day to estimate erg . We then examine the evolution of 
cluster properties to constrain the value of Qq- Specifically, 
we present predictions for the mass function, the temper- 
ature function and the distribution of S-Z decrements at 
z ~ 0.3 — 0.5. This range of redshifts is a rich area of obser- 
vational work. Ongoing programmes include mass measure- 
ments from weak gravitational lensing using HST and large 
format CCDs (Ellis, private communication; Kaiser, private 
communication); estimates of the cluster temperature func- 
tion (Henry, private communication); and S-Z source counts 
and redshift distributions from ground-based radio tele- 
scopes and, if approved, from the COBRA-SAMBA satellite 
mission. 

Related calculations have been carried out by Evrard 
(1989), Frenk et al. (1990), Bond & Myers (1992), Lilje 
(1992), Oukbir & Blanchard (1992), Bahcall & Cen (1993), 
Hanami (1993), White et al. (1993), Bartlett & Silk 
(1993,1994), Barbosa et al. (1995), Hattori & Matsuzawa 

(1995) , Viana & Liddle (1995), amongst others. Our work 
differs from these earlier studies in various respects. Firstly, 
we consider a wider range of cosmological models than most 
previous analyses, particularly models with low values of Jlo, 
with and without a cosmological constant. Secondly, we use 
some of the largest N-body simulations ever carried out in 
order to verify the accuracy of our analytical calculations. 
Finally, we place special emphasis on making detailed pre- 
dictions in a manner that can be compared to observational 
data with a minimum of assumptions. Since high-redshift 
data are not yet sufficiently complete, we do not attempt to 
estimate fio in this paper, but present instead a number of 
distributions that will enable this estimate to be made as 
the data become available. 

This paper is structured as follows. Section |2| reviews 
the Press-Schechter model for the abundance of clusters, ex- 
tended to open and flat cosmological models {Qo < 1 with 
Ao = or Ao = 1 — Slo). The predictions of this model are 
compared with a set of 192^ particle N-body simulations in 
Section |^. In Section ^ we review the existing cluster abun- 
dance data and normalise the analytical models using the 
cumulative cluster temperature function estimated from the 
data compiled by Henry & Arnaud (1991). The evolution of 
the cluster mass and temperature functions, the S-Z source 
counts and the dependence of t hese quantities on fio and Ao 
are predicted in Sections - 5^. Finally, we discuss the 
importance and implications of these results in Section ^. 



2 MODEL 

An analytical expression for the comoving number density of 
dark matter halos of mass M in the interval dM, originally 
derived by Press & Schechter (1974) (see also Bond et al. 
1991), is 



© 0000 RAS, MNRAS 000, 000-000 



Cluster Evolution 3 



1.7 



1.65 



l.i 



150 - 



<i 100 - 



50 



A-0 



— n+A=i 



A-O 
Q I A=l 



0.2 



0.4 0.6 

n 



0.8 



Figure 1. Upper panel: critical tlireshold for collapse, Sc, as a 
function of Q, in tlic splierical collapse model. Results are plotted 
for open models with A = (dotted line) and fiat models with f!+ 



A = 1 (dashed lines). Lower panel: the virial density of collapsed 
objects m units ot tne critical aensity. ine aottea ana aasnca 
lines are as in the upper panel. 
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(2.1) 



where p is the present mean density of the universe, and 
cr(M) the present, linear theory rms density fluctuation in 
spheres containing a mean mass M . The evolution with red- 
shift is controlled by the redshift dependent density thresh- 
old 5c (2). For the case of fJo = 1, the conventional choice of 
this threshold is (5c(2) = 1.686(1 + z). This is the value, ex- 
trapolated to the present using linear theory, of the overden- 
sity of a uniform spherical overdense region at the point at 
which the exact non-linear model predicts that it should col- 
lapse to a singularity. This threshold, along with the choice 
of top-hat filtering to define (j(M), gives a mass function 
that agrees remarkably well with the results of N-body sim- 
ulations (e.g. Efstathiou et al. 1985, Lacey & Cole 1994). For 
f2 7^ 1 one can use the spherical collapse model to derive a 



general expression for the density threshold. Expressing the 
threshold as = <5c(0)/D(2;,f2o, Ao) where D(z, r2o , Ap 



is the line ar growth factor normalised to unity at z = (Pee- 



bles 198C), we find that 5c (0) has only a weak dependence 



on for both open models with A = and flat models with 
n -f- A = 1 (see Fig |l]). The details of the A = calculation 
can be found in Lacey & Cole (1993; see also Maoz 1990), 
and the -I- A = 1 result is derived here in the Appendix 
(see also Lilje 1992, White et al. 1993 and Kochaneck 1995). 

^In order to convert the mass function obtained from 

( ^.l| ) to a temperature function we assume that the gas is 
isothermal. In this case. 
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(2.2) 



where Ac is the ratio of the mean halo density, to the criti- 
cal density at that redshift, (5 is the ratio of specific galaxy 
kinetic energy to specific gas thermal energy and X is the hy- 
drogen mass fraction which we take to be X = 0.76. Recent 
work by Navarro et al. (1995) shows that equation (2.2) is ac- 
curately obeyed in N-body/hydrodynamic simulations of the 
formation of clusters in universes with f2o = 1- Their simu- 
lations predict an X-ray luminosity-weighted /3 = 1.07±0.05 
for individual clusters, in fair agreement with observational 
determinations (e.g. Forman & Jones 1990). We adopt the 
value /? = 1 throughout this paper, but our results may 
be modified for different choices simply by rescaling all pre- 
dicted temperatures by /3~^. 

The density contrast Ac is computed from the spher- 
ical collapse model assuming that the cluster virialises at 
the redshift at which we view it. Its dependence on f2 (and 
therefore on redshift) is given in the lower panel of Fig |l[ 
The assumption that clusters form at the redshift at which 
we view them is a good approximation in the fio = 1 model 
since in this case halos are continuously accreting material. 
The value Ac = 178, appropriate to = 1, has recently 
been shown to separate well the interior of the virialised 



halo from the surrounding infalling material ( dole & Lacey 



1995). In low-f2 models the accretion rate onto a halo de- 
clines after z ~ 1/(0,^^ — 1) and so little evolution in the 
density and temperature of the bulk of the gas takes place 
at low redshift. The density in the simple spherical collapse 
model, Pcrit(l + zf{Q,o/VL{z))l^c (where p% it is the critical 
density at the present day), does not accurately reproduce 
this behaviour but the high resolution simulations needed 
to determine it have yet to be carried out. Nevertheless, for 
the range of redshifts and f2os that we consider, the simple 
model is a good approximatioii; 



The Sunyaev-Zel'dovich (Sunyaev & Zel'dovich 1972 



effect is produced by the inverse Compton scattering of 
cosmic microwave background (CMB) photons off high en- 
ergy electrons in the intracluster gas. This process dis- 
torts the blackbody spectrum by shifting microwave pho- 
tons to higher energies. At long wavelengths (longward of 
A = 1.37 mm for Tcmb = 2.726K), the cluster produces a 
negative fiuctuation in the surface brightness of the CMB, 
while at shorter wavelengths it produces a positive fiuctu- 
ation. In the long wavelength regime the microwave back- 
ground decrement from the cluster is given by 
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AT 



-2y, 



(2.3) 



where y is the integral of the electron pressure along a line- 
of-sight through the cluster, 



y 



ncfJT 



(^) 



dl. 



(2.4) 



Here rie is the number density of electrons, rrie the electron 
mass, and (tt the Thompson cross-section. We can define an 
effective angular cross-section, Y , by integrating y over the 
projected area of the cluster and dividing by the square of 
the angular diameter distance, r^, 



Y 



y dA. 



(2.5) 



This quantity has a simple physical interpretation. In the 
long wavelength regime 2Y is simply the efi'ective angular 
area of the microwave background obscured by the clus- 
ter. An unresolved observation of a cluster with a radio 
telescope of effective beam area, Abeam, would, in the long 
wavelength regime, measure a signal AT/T = — 2y/Abcam- 
If the cluster is resolved, then the signal depends on the 
density profile of the cluster. Assuming a surface den- 
sity proportional to , as in an isothermal sphere, then 
AT/r = -2y/(Acius Abeam)^'^^ whcrc Adus IS the angular 
cross-section of the virialised cluster. Using the virial radius 
to define the edge of the cluster, we can then write 



= 64.1 



/.(5X-,3)(^)( 



l + z) 



V178/ VlOOh-iMpc 



(2.6) 



The effective angular area, Y, is also easily related 
to the CMB flux absorbed (or re-emitted) by the clus- 
ter at other frequencies. Defining a dimensionless frequency 
X — hpu/kTcMB = Ao/A, where hp is Planck's constari t and 
Ao = 5.28mm for Tcmb = 2.726K ( [Mather et al. 1994[ ), the 
emitted flux is 



S.(x) = S'i^^^ {x) Q{x) Y 



= 2.29 X 10* 



Y 



mjy, 



where 



Q{x) = 



(2.7) 
(2.8) 

(2.9) 



tanh(x/2) 

^m~^ Hz~^. In the long wavelength 



and 1 mJy= IQ-^^ Js 
limit, X ^ Q, Q{x) —2 and this reduces to the result 
stated earlier, AT/T = S^/S^^^ = -2Y. At x ^ 3.83 (A = 
1.37mm) Q{x) = and at higher frequencies the cluster 
appears as a source of emission 



The effective cross-section, Y, deflned by equations (2.4) 
and (2.5), is proportional to the total mass and average tem- 
perature of the intracluster gas. Speciflcally, for an isother- 
mal intracluster gas. 
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X 10-'h/icM(l + X) 
'kT 
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IQiSh-iM 



\Tev) 



rd 



lOOh-iMpc 



(2.10) 



where /icm is the fraction of the cluster mass represented by 
the hot intracluster gas. We adopt /3 = 1 and /icm = 0.1 but 
our results can readily be rescaled to other choices. Note that 
here it is appropriate to use the mass-weighted temperature 
for which Navarro et al. (1995) flnd /3 — 1.2 for their flo — 1 
simulations. 



3 COMPARISON WITH N-BODY 
SIMULATIONS 

In order to assess the accuracy of the Press-Schechter mass 
distribution in the regime of interest - the mass scale of rich 
galaxy clusters - we compare the model predictions with 
the abundance of dark matter clumps found in a new set 
of large cosmological N-body simulations (Cole, Frenk & 
Weinberg in preparation). The sim ulation s were performed 



with the AP'^'M code of Couchman (1991) using 192'=' « 7 x 



10 particles in a periodic box of size /box = 345.6h-iMpc. 
The particle mass was Mp = 1.64 X 10^^noh"^MQ and the 
force resolution e — 180h~^kpc, where e is the equivalent 
Plummer potential softening parameter. Two sequences of 
simulations were carried out with different values of 57o, a 
sequence of open models with Ao = and a sequence of 
flat models with flo + Ao = 1. Both were normalised to 
have as = 0.55 f^o^"'* so as to reproduce, approximately, the 
observed abundance of galaxy clusters (White et al. 1993). 
The number of timesteps required to evolve the simulations 
accurately from the linear regime to the present epoch was 
approximately 100/Slo- The initial linear power spectrum 
was the same in all simulations, a CDM-like spectrum with 
scale parameter F — f2oh = 0.25. This value of F is suggested 
by observations of large-scale galaxy clustering (eg. Maddox 
et al. 1990). Here we will consider only three representative 
models, one with fio ~ 1, another with Qo = 0.3 and Aq — 0, 
and a third one with flo = 0.3 and Aq = 0.7. 

We identifled groups of particles in the simulations us- 
ing two different algorithms. The first was the standard 
friends-of-friends algorithm (Davis et al. 1985) with link- 
ing length bi times the mean interparticle separa tion; t he 



second wa s the spherical overdensity algorithm (Lacey & 



Cole 1994) with density contrast ftp. We chose values of 



characteristic of virialised objects predicted by the spheri- 
cal collapse model. For the friends-of-friends linking length 
we took bi — 0.2 in the 57o = 1 simulation and scaled as 
bi oc (Ac/0,)~^^''^ for the other cases. In both the simula- 
tions and the Press-Schechter model we converted mass to 



temperature using relation (2.2). 

Fig 1^ compares cluster abundances as a function of kT 
at z = and z — 0.5 in the N-body simulations and in 
the Press-Schechter model. Overall, the analytical predic- 
tions are in excellent agreement with the N-body results 
and reproduce the near exponential fall-off of the temper- 
ature functions very accurately. The threshold density, Sc, 
has not been treated as a free parameter, but has instead 
been fixed at the value prescribed by the spherical collapse 
model. For the f2o = 1 case at z — this value is es- 
sentially perfect and produces a temperature function mid- 
way between that obtained using the friends-of-friends and 
spherical overdensity algorithms. For z = 0.5, the spherical 
overdensity algorithm finds somewhat fewer high tempera- 
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Figure 2. Comparison of the temperature functions predicted by the Press-Schechter distribution and the results of N-body simulations. 
For each of the three cosmological models (with parameters given at the top of each panel), temperature functions are plotted for z = 
and z = 0.5. The Press-Schechter predictions, which are normalised according to the spherical collapse model, are shown by solid lines. 
The simulation results are plotted with dashed lines for clusters identified with a friends-of-friends algorithm and with dotted lines for 
clusters found with the spherical overdensity algorithm. Over the full range of abundances and redshifts shown the agreement between 
theory and simulations is very good. 



ture clusters than either the friends-of-friends algorithm or 
the Press-Schechter model. This may partially reflect the 
relatively low resolution of the N-body simulations as the 
most massive groups contain only about 100 particles in the 
f2o = 1 simulation at z = 0.5. In the two Qo = 0.3 cases, 
the Press-Schechter model matches the very slow evolution 
of the halo abundances in the simulations extremely well. If 
one were to adjust 5c upwards by approximately 4%, then 
the N-body results would be reproduced almost perfectly 
over the temperature range plotted in Fig ^. However, such 
a small adjustment is negligible compared to the uncertainty 
in the abundance of real clusters. We note that for the more 
common, lower temperature clusters not shown in this fig- 
ure, the Press-Schechter distribution predicts a significantly 
higher abundance than found in the simulations. We con- 
clude that, for the range of masses of interest, the agree- 
ment with the numerical results is sufficiently good that we 
can confidently use the Press-Schechter model for detailed 
calculations. 



4 NORMALISATION OF THE MODELS 

We now compare the predicted and observed distributions 
of cluster X-ray temperatures in the local universe in or- 
der to obtain an estimate of erg, the rms mass fluctuation 
in spheres of radius 8h~^Mpc. The observational data we 
use is the complete flux-limited sample of 25 clusters com- 



piled by Henry & Arnaud (1991). We compare our estimate 



with earlier determinations of ag from both X-ray and op- 
tical data (Henry & Arnaud 1991; White et al. 1993; Viana 
& Liddle 1995) and discuss the reasons why these earlier 
determinations gave slightly different values of erg. 

The stepped curve in Fig ^ shows the following sim- 
ple estimate of the cumulative cluster temperature function 
derived from the Henry & Arnaud data. 



iV(> kT) = 



kTi>kT 



1/Vm 



(4.1) 



where Vmax,i is the maximum volume in which the cluster 
could be detected given the flux limit and geometric bound- 
aries of the survey. We have chosen to present the results in 
cumulative form in order to avoid binning the data. Since the 
cluster abundance falls very rapidly with increasing temper- 
ature, the differential temperature function averaged over 
each bin can be significantly larger than the underlying un- 
binned distribution. Henry & Arnaud (1991) estimated the 
differential temperature function by weighting the number 
of clusters in each temperature bin with the average value of 
Vinax.i in that bin. This estimator is only equivalent to ( |4.l| ) if 
Vinax.i is the Same for all of the clusters in each bin. However, 
since the bins have non-zero width and the Lx-Ix relation 
has considerable scatter, there is also considerable scatter 
in the individual Vinax.i values in each temperature bin. For 
the Henry & Arnaud dataset this is a large effect and would 
have led them to underestimate the cluster abundance by 
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Figure 3. Predicted and observed X-ray temperature functions. The left-iiand panel shows model predictions for T = 0.25 and f2o = 1 
(smooth solid line), Oq = 0.3, Aq = (dotted line) and Oq = 0.3, Aq = 0.7 (dashed line). The models are normalised by fitting to the 
observed data as described in the text. The observed temperature function, shown as a solid line with steps, was derived from the data 
compiled by Henry & Arnaud (1991). The error bars at the three temperatures where the models were fit were obtained from a bootstrap 
analysis. The best-fit values of a-g are given in the figure. The long-dashed line shows the cumulative temperature function obtained by 
Edge et al (1990). The hatched box represents the range of normalisations derived by White, Efstathiou & Frcnk (1993) from a similar 
theoretical analysis of combined X-ray and optical data. The right-hand panel gives model predictions for f2 = 1, erg = 0.5 and three 
values of F: 0.15 (solid line), 0.25 (dotted line), and 0.50 (dashed line). The observational data and the hatched box are the same as in 
the left-hand panel. 



approximately a factor of 4. However, (Henry private com- 
munication) in computing each Knax.i, a spurious factor of 
4.2 entered their calculation and this largely compensated 
for the bias in the estimator. Thus the published Henry & 
Arnaud (1991) temperature function is in reasonably good 
agreement with our estimate. It is also in reasonable agree- 
me nt wit h the cumulative temperature function of Edge et 
al. (199C), reproduced as the long dashed line in Fig ^. 

To estimate ag, we fit our model predictions to our 
estimate of the temperature function. First, we calculate 
the statistical errors in the temperature function through a 
boostrap procedure. At three selected temperatures, Ti, we 
computed (using equation 4.1) the temperature function for 
a large number of samples of 25 clusters constructed by ran- 
domly selecting, with replacement, from the original list of 
25 clusters in Henry & Arnaud's compilation. The error bars 
in Fig. ^ show the resulting 1 — a ranges in the bootstrap 
distribution of logjQ A'^boot (> kT), at our three chosen tem- 
peratures. This distribution was also used to compute the 



covariance between the estimates at different temperatures. 
Averaging over the bootstrap samples we thus obtain the 
covariance matrix, 

Crj = (e.e,), (4.2) 

where = log^p A'boot(> fcTi)-log^o A'^data(> kTi). We then 
fit the model temperature functions by minimising 



X 



(4.3) 



where = log^g A'niodei(> kTi) - logig A^data(> kT^). If 
the data points are uncorrelated dj is diagonal and equa- 
tion (4.3) reduces to the normal definition of x^- We find 
that there are significant correlations between the esti- 
mates of logjQ N{> kT) at the three selected temperatures, 
but that the models which minimise ^^re insensitive to 
whether the correlations are treated as above or simply ig- 
nored. 

The best-fit models are shown in Fig ^ for fio = 1 and 
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Slo = 0.3, with and without a cosmological constant. The 
models in the left-hand panel all have the same CDM-like 
power spectrum with F = 0.25. Because of the very sensi- 
tive dependence of cluster abundance on spectrum normali- 
sation, erg can be estimated with high precision even though 
the errors in the empirical cluster abundance are quite large. 
The formal error on as from the fits in Fig 3 is ±2%, but 
this is likely to be an underestimate of the true uncertainty 
because, with only 25 clusters, the errors in the tempera- 
ture function are unlikely to be Gaussian and, in addition, 
systematic errors are likely to be significant. 

The right-hand panel of Fig 3 shows the effect of allow- 
ing r to vary for the case of fio = 1. Values of F < 0.25 
produce more large scale power and more very hot clus- 
ters. This produces a temperature function with a shallower 
slope, which is in better agreement with the observations. 
(The best fitting model has F — 0.03 and ug = 0.51, but this 
is not a reliable way of constraining F because the high and 
low temperature points are anticorrelated, and this creates a 
large uncertainty in the slope of the temperature function.) 
Fortunately, the best fitting values of erg are quite insensitive 
to the adopted value of F, because as F varies it is the slope 
of the temperature function that varies with the pivot point 
remaining in the middle of the range of rich cluster temper- 
atures. This simply reflects the fact noted above that the 
average mass within a sphere of radius 8h~^Mpc in an un- 
perturbed universe with f2o ~ 1 is very close to the mass of 
a rich galaxy cluster. 

In summary, we flnd that the observed temperature 
function is well fit by all our CDM-like models if 



erg = (0.50 ± 0.01)^0 ° *^+°-'''"'' for Ao = 



(4.4) 



and 

erg = (0.50 ±0.01)^0 



0. 53-1-0. ISflo 



for fio + Ao = 1. (4.5) 



The quoted statistical uncertainty corresponds to the stan- 
dard deviation of the as values obtained from fitting the 
individual bootstrap catalogues. A more realistic estimate 
of the error would take into account the variation of erg due 
to the uncertainty in the measured X-ray temperatures and 
the values of both 5c and /3. The error in 5c we estimate from 
Fig. 2 to be < 4%. The uncertainty in the appropriate /? we 
take to be 10% (see Navarro et al. 1995). Note that adopt- 
ing /3 = 1.1 leads to values of erg which are 6% larger than 
those quoted above. Combining these sources of error we es- 
timate the overall uncertainty in erg to be approximately 8%, 
which is four times as large as the statistical error shown in 
equations (4.4) and (4.5). 

The values of erg that we infer from the X-ray data are 
systematically lower than those obtained by White et al. 
(1993) who found erg = 0.57^^°-^^. Our modelling of the 
cluster abundances is almost identical and the difference in 
our inferred values of as results almost entirely from the 
different observational data that we fitted. White et al. esti- 
mated the mass of clusters with abundance 4 x 10~®h'^Mpc~'^ 
in two different ways. The first and larger estimate was based 
on the median velocity dispersion of rich Abell clusters. The 
second was based o n the cumulative temp eratu re functions 
of Henry & Arnaud ( [L991| ) and Edge et al. (|l990| ). The range 
spanned by these two estimates, expressed as a temperature 
rather than a mass using equation (2.2), is indicated by the 
hatched box in Fig ^ The lower estimate of the X-ray tem- 



perature is slightly higher than our analysis of the Henry 
& Arnaud data implies. This difference is due to the fact 
that the White et al. estimate was inferred indirectly from 
the differential temperature function presented by Henry & 
Arnaud, not directly from the cumulative temperature func- 
tion. (If we adopt as = 0.570,^'^'^^ then our model tempera- 
ture curves pass through the shaded box in agreement with 
White et al. ). 



Henry & Arnaud (1991) found erg = 0.59 ± 0.02 for 



fio = 1 and a power-law fluctuation spectrum. They used 
essentially t he sa me data as us, but they adopted f5 — 1.2 
in equation (2.2), while we have assumed /3 = 1. Thus, for 
the same value of erg, their model temperature functions are 
shifted to lower temperatures by ~ 0.08 in logj^Q(fcr), and 
this largely accounts for the difference in the inferred values 
of erg. Viana & Liddle (1995) obtained as = 0.6 for ilo = 1, 
with a dependence on flo close to that which we find, by 
fitting only to the abundance of clusters at 7keV. It may 
be seen from Fig ^ that fitting the cumulative temperature 
function just at this temperature yields a higher value of erg 
than the one obtained from fitting all the X-ray data. Given 
the small number of clusters with temperatures as great as 
7keV, fitting over a wider temperature range seems more 
appropriate. 

In summary, our modelling of the cluster temperature 
function is co nsiste nt with those of Henry fc Ar naud (1991), 
White et al. ( |l993| ) and Viana & Liddle ( |l995| ). The range 
in the values of erg deduced by these authors arises from the 
different data poi nts th ey chose to fit and, in the case of 
Henry & Arnaud (1991), from the value of /? they used to 
relate virial to gas temperature. Assuming /? = 1, we con- 
clude that the X-ra y da ta are best fit by the values of erg 
given by equations (4.4) and (1.5). Since the abundance of 
Abell clusters of richness class > 1 is 8 x 10~®h'^Mpc~'^, our 
results imply that the median 1-D virial velocity dispersion 
of these clusters should be approximately 650km s^^. This is 
smaller than the median values of around 800km s~^, found 
in the compilations of Zabludoff, Huchra & Geller (1990) 
and Girardi et al. (1993). These differences may be under- 
stood if, as argued by Frenk et al. (1990) and others, the 
higher velocity dispersion estimates are artificially boosted 
by contamination from infalling groups around the cluster. 
Alternatively, these larger velocity dispersions (and a higher 
value of erg) would be consistent with the X-ray data if the 
intracluster gas were significantly cooler than the virial tem- 
perature, but this would require /3 > 1.5. Such large values 
of /3 are not supported either by the data or by recent hy- 
drodynamic simulations of cluster formation (Navarro et al. 
1995, Evrard et al. 1995). 



5 CLUSTER EVOLUTION 

Having fixed the normalisation of the models by requiring 
that they should match the local abundance of rich clus- 
ters, we now consider their evolutionary properties. Specifi- 
cally, we calculate the redshift dependence of the mass func- 
tion. X-ray temperature function and the distribution of S- 
Z decrements. For each model, the evolution of the clus- 
ter mass function (equation ( [2.l| ) is governed by the linear 
growth factor, D{z, fio, Ao). The evolution of the X-ray tem- 
perature and S-Z effect depends, in addition, on the evolu- 
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Figure 4. Predicted evolution of the cluster mass function. The 
comoving number density of clusters per (h~^Mpc)'^ with masses 
larger than M is shown as a function of Af . Solid lines correspond 
to n = 1; dotted lines to an open model with Qq = 0.3; and 
the dashed lines to a flat model with = 0.3 and Aq = 0.7. 
Predictions for z = 0, z = 0.33 and z = 0.5 are plotted. There 
is relatively little evolution in the Qg < 1 cosmologies but, in an 
Qq = 1 model, the abundance of clusters declines precipitously 
with redshift. 
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tion of the virial density Aciio/n{z), which determines the 
relation (^^) between mass and temperature. 



5.1 Evolution of the cluster mass function 

The growth of fluctuations continues at a rapid rate at re- 
cent times if f2o = 1 and very little if Qq = 0.3. As a result, 
the rich cluster mass function evolves dramatically between 
z = 0.5 and z = if S7o = 1, but much less so if flo is low, 
whether or not the cosmological constant is zero. This evolu- 
tion is illustrated in Fig ^. For f2o = 1, the comoving number 
density of clusters of virial mass M = 3.5 X lO^h-^Mo, typ- 
ical of Abell clusters of richness class R > 1, declines by a 
factor of 30 between z — and z — 0.33 and, by 2; = 0.5, it 
is tiny. By contrast, for Qq — 0.3, Ao — 0.7, the abundance 
of clusters of this mass has only dropped by a factor of ~ 7 
below the present day abundance even at z = 0.5; if Ao = 
the decline is even slower. 

The strong f2o-dependence of the rate at which the mass 
function evolves is reflected in the expected redshift distribu- 
tions of massive clusters (A/ > 3.5 x lO^*h~^M0), illustrated 
in the top panel of Fig ^. For f2o = 1 the distribution peaks 
sharply at very low redshift, whereas for low flo a broader 
peak is displaced to higher redshift. The effect of Ao is to 
move the peak back to a somewhat lower redshift, reflecting 
the slightly later epoch at which structure ceases to grow in 
non-zero Ao cosmologies. 



Figure 5. Upper panel: Redshift distribution of massive clusters 
{M > 3.5 X lO^''h~^M0) in different cosmological models. The 
ordinate gives the probability distribution of clusters per unit 
redshift interval. Lower panel: Number counts of clusters with 
mass M > 3.5 X lO^*h~^M0 out to a given redshift. The ordinate 
gives the count per unit area on the sky. In both panels solid lines 
correspond to CIq = 1, dotted lines to = 0.3 and dashed lines 
to Oq = 0.3, Ao = 0.7. The models are normalised by the value 
of as for which the predicted temperature function at 2 = best 
fits the data. The low-Ho cosmologies produce significantly more 
clusters at high redshifts than the Qq = 1 model. 



Integrating over the redshift distributions yields the 
number count of clusters per unit area on the sky. The 
change in the volume element corresponding to a fixed red- 
shift interval enhances the differences between the high and 
low-f2 models. At redshift z = 0.33, the volumes per unit 
redshift are in the ratios 1 : 1.23 : 1.7 for f^o = 1, Jlo = 0.3, 
Ao = and Qq = 0.3, Ao = 0.7 respectively. At z = 0.5 
the corresponding ratios increase to 1 : 1.35 : 2.02. Thus, for 
Slo = 1, we expect to find only 0.015 clusters per square 
degree with mass greater than 3.5 x lO'^*h~^M0 out to 
z = 0.5 and virtually none at higher redshifts. By contrast, 
for 0,0 — 0.3, we expect to find more than 10 times as many 
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clusters above this mass with z < 0.5. Note that predictions 
for low-f2 models are relatively insensitive to the value of 
Ao, with only about 50% more clusters predicted to exist in 
the non-zero Ao model. This factor results largely from the 
difference in the volume elements. 

The cluster mass that enters into equation (2.1) and 
Figs ^ and ^ is the virial mass, ie the mass contained within 
a sphere of mean overdensity Ac. In practice, gravitational 
lensing measurements give the mean projected surface den- 
sity, E(7i), within a radius, R, which is typically less than 
lh~^Mpc. To relate these two masses requires a model of 
the cluster mass distribution. Recent high resolution N-body 
simulations show that, over most of the cluster, the dark 
matter density profile is well approximated by an isothermal 
profile (Navarro et al. 1995; Cole & Lacey 1996). Although 
in principle it would be straightforward to use the analytical 
fit to the N-body profile given by Navarro et al. (1995), for 
our purposes the isothermal profile approximation is quite 
adequate. In this case, mean surface density is related to 
virial mass by 



Mvir = 



4 / GR^T.{Rf 



H 



1/2 



(5.1) 



where R is any radius inside the virial radius, H is Hubble's 
constant at the redshift of the cluster and the overdensity 
Ac is given in Fig 1 as a function of fi. This formula may 
be used to relate the surface mass density estimated from 
weak gravitational lensing analyses to the virial mass used 
throughout this paper. Note, however, that we have implic- 
itly assumed that any foreground or background mass makes 
a negligible contribution to the lensing signal. 



5.2 Evolution of the X-ray temperature function 

The evolution of the X-ray temperature function depends 
both on the growth factor, 73(2, f2o, Ao), and on the virial 
density, /S.c^o/^{z). Thus, at high z, the abundance of clus- 
ters of a given temperature is determined by the balance 
between the overall decline in the population of virialized 
clusters and the lower mass associated with each tempera- 
ture. If f2o = 1, the first factor is dominant and the tem- 
perature function declines precipitously (see Fig ^) . By con- 
trast, in an open, fio = 0.3, universe, the two effects nearly 
cancel out and the temperature function remains virtually 
unchanged at least out to 2; = 0.5. A flat, Q,o = 0.3, universe 
is intermediate and, in this case, the temperature function 
declines slowly with redshift. 

The redshift distribution of clusters hotter than 5 keV 
is shown in Fig |^. As was the case for the mass function, the 
f2o = 1 model produces large clusters predominantly at low 
redshifts. The Q,o = 0.3 models, on the other hand, give rise 
to extended redshift distributions. Again, the lower redshift 
at which structure "freezes-out" when a Ao term is included 
produces somewhat stronger evolution in this case compared 
to an open cosmology. When integrating over redshift, the 
effect of different evolutionary rates is enhanced by the larger 
volume in low density universes. As a result, the number 
counts per unit area of sky, displayed in the lower panel of 
Fig ^ depend strongly on the value of f2o. For our choice 
of parameters, the counts in the two fio = 0.3 models are 
very similar and the total number of clusters hotter than 




0.6 0.8 
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Figure 6. Predicted evolution of the cluster X-ray temper- 
ature function. The comoving number density of clusters per 
(h~^Mpc)^ hotter than kT is shown as a function of kT. Solid 
lines correspond to f2 = 1; dotted lines to an open model with 
Co = 0.3; and the dashed lines to a flat model with Qo = 0.3 
and Ao = 0.7. Predictions for z = 0, z = 0.33 and z = 0.5 are 
plotted. There is little evolution in the Co < 1 cosmologies but in 
an Co = 1 model, the abundance of clusters declines precipitously 
with redshift. 



5 keV at redshifts less than 0.5 is about 9 times higher in 
these models than in the Q.o = 1 case. The expected number 
counts for different values of f2o are shown in Fig ^ for both 
open and fiat models. 

We emphasize that the results shown in Figs |^ - ^ are 
very sensitive to the normalisation of the fluctuation ampli- 
tude. Our adopted values of erg were fixed by requiring that 
each model should agree well with the present day X-ray 
temperature function. With this particular choice, the in- 
clusion of a Ao term turns out to make very little difference 
to the predicted abundance of hot clusters out to redshift 
of 0.5. However, the number counts out to this redshift do 
discriminate well between different values of JIq. 



5.3 Evolution of Sunyaev-Zel'dovich effect 

As discussed in Section 2, we characterize t he S- Z effect by 
the y-function given in equations (2.5) and (2.10). (Here, we 
take /icM = 0.1 and /9 = 1, but our results are readily scaled 
to other values using Y oc /icm/3~^.) The evolution of the Y- 
function is similar to that of the temperature function, but 
since Y oc M^/^ whilst T oc M^^^ , the detailed behaviour is 
slightly different. More importantly, the appearance of the 
diameter-distance in equation (2.5) causes Y to drop off with 
redshift more rapidly than kT. 

Figure ^ confirms the rapid decline in the y-function 
with redshift, resulting from the r^'^ term, in all our cosmo- 
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Figure 7. Upper panel: Redshift distribution of hot clusters 
(fcT > 5keV) in different cosmological models. The ordinate gives 
the probability distribution of clusters per unit redshift interval. 
Lower panel: Number counts of clusters hotter with kT > 5keV 
out to a given redshift. The ordinate gives the count per unit area 
on the sky. In both panels solid lines correspond to fip = 1, dot- 
ted lines to Qq = 0.3 and dashed lines to Qq = 0.3, Aq = 0.7. The 
models are normalised by the value of eg for which the predicted 
temperature function at z = best fits the data. The low-f2o 
cosmologies produce significantly more clusters at high redshifts 
than the Oq = 1 model. 

logical models (which, as before, are normalised using the lo- 
cal cluster X-ray temperature function). Evolution is in the 
same sense as evolution in the cluster mass and X-ray tem- 
perature functions. Also, for z > 0, the differences between 
the y-functions of our various models are of similar magni- 
tudes. To understand why this is so we need to consider the 
two factors mentioned above. In the Qo — 1 model, in which 
the fluctuation growth factor is still changing rapidly at low 
redshifts, M^^^ decreases faster with redshift than either M, 
which is relevant to the mass function, or M^^'^, which is 
relevant to the temperature function. Thus, for Qq = 1, the 
y-function evolves more rapidly than the mass or temper- 




Figure 8. Predicted number counts of clusters hotter than 5 keV 
out to z = 0.33 and 0.5, as a function of Qq. The solid and dotted 
lines correspond to open models (Aq = 0) and the dashed lines to 
flat models (Aq = 1 — Oq). The abundances are sensitive to the 
choice of trg, but the ratio of the numbers expected in different 
cosmologies is less so. 



aturc functions. This is partly offset, in the low-Slo models, 
by their larger angular-diameter distance out to a particular 
redshift. 

Because of the strong evolution apparent in Fig ^, the 
redshift distribution of "S-Z bright" clusters is highly peaked 
at low redshift in all our models. This distribution is plotted 
in Fig [lo| for clusters with Y > 10^'' h arcmin'^. The corre- 
sponding number counts as a function of Qq are shown in 
Fig ^ The counts are dominated by local clusters and, since 
all the models are normalised to the local temperature func- 
tion, the counts are similar in all cases. Nevertheless, ~ 60% 
more Y > 10~'^h arcmin^ clusters are expected if f2o ~ 0.3 
than if fio = 1- The splitting between the various cosmolo- 
gies increases if the Y threshold is reduced, as illustrated in 
Fig As the threshold is lowered, the counts probe higher 
redshifts where the predictions are increasingly sensitive to 
Qq. Thus, if the threshold is taken to be 10~''h arcmin^, 
over twice as many clusters are expected per unit area if 
Qo = 0.3 than if fio = 1- As with the counts as a function 
of mass or X-ray temperature, the counts as a function of 
Y are insensitive to Aq. The redshift distribution of a sam- 
ple of Y-selected clusters gives a more promising method of 
constraining Q,o- For clusters with Y > 10~"^h arcmin^ the 
median redshift (see Fig. 10) in the two ilo = 0.3 models is 
approximately z = 0.136, which is a little over twice that 
of the r^o = 1 model. Thus optical follow up on a relatively 
small sample of Y-selected clusters could easily distinguish 
between these two models. 
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Figure 9. Predicted evolution of tlie cumulative cluster S-Z Y- 
function. The ordinate gives the comoving number density of clus- 
ters per (h~^Mpc)''. The solid lines are for Qg = 1; the dotted 
lines for fig = 0.3 and the dashed lines for Qq = 0.3 and Aq = 0.7. 
Predictions for z = 0.05, z = 0.33 and z = 0.5 are plotted. Yhe 
y-function evolves significantly in the three cosmologies, but the 
evolution is strongest for Qq = 1. 



6 DISCUSSION AND CONCLUSIONS 

We have used the Press- Schecther formalism to investigate 
the evolution of the population of rich galaxy clusters. Our 
work extends and complements earlier work using this for- 
malism by Evrard (1989), Henry & Arnaud (1991), Lilje 
(1992), Oukbir & Blanchard (1992), Hanami (1993), White 
et al. (1993), Barbosa et al. (1995), Hattori & Matsuzawa 
(1995), and Viana & Liddle (1995), amongst others. We first 
verified that the Press- Schechter formula predicts the correct 
abundance of rich clusters at 2; = and 2: = 0.5 by compar- 
ing with the results of large cosmological N-body simula- 
tions. The agreement is excellent, at least for the models we 
have considered: CDM-like cosmologies with spectral shape 
parameter, F = 0.25 and with (a) Qo = 1, (b) Slo = 0.3, 
Ao = and (c) ilo = 0.3, Ao = 0.7. 

We reconsidered the problem of normalising the ampli- 
tude of mass fiuctuations on cluster scales by reference to 
the present day abundance of rich clusters. From our own 
rederivation of the X-ray temperature distribution of clus- 
ters, using Henry & Arnaud's (1991) data, we found the 
following values for the rms density fluctuation in spheres of 
radius 8h~^Mpc, ag: 



(78 = (0.50 ± 0.04:)Qo°'^''+°-^°"° for Ao = 



(6.1) 



and 

as = (0.50 ± 0.04)^0 



0.53-1-0.1300 



for no + Ao = 1. (6.2) 



Figure 10. Redshift distribution of S-Z bright clusters (Y > 
10~^h arcmin^) in different cosmological models. The ordinate 
gives the probability distribution of clusters per unit redshift 
interval. The solid line corresponds to f2o = 1, dotted line to 
f2o = 0.3 and dashed line to Hg = 0.3 and Ao = 0.7. 



the power spectrum and, for other values of Qo, there is only 
a very weak dependence. Note that erg is only slightly larger 
if the A term is non-zero. For fio > 0.2, the difference be- 
tween the flat and open models is always less than 10%. Note 
also that since th e rms fluctuation in the bright galaxy dis - 
tribution is 0.96 (Maddox, Efstathiou, & Sutherland 199C) 
in spheres of radius 8h^ Mpc, our estimates of as imply that 
the biasing parameter deflned as fo = 0.96/a"8 is greater than 
unity for open models with Qq > 0.23 and for flat models 
with Qo > 0.27. Models with fio smaller than this require 
antibiasing, that is they require bright galaxies to be less 
clustered than the mass. Our estimates of as differ slightly 
from those obtained by previous authors and, in Section 4, 

we discussed in detail the reasons for these differences. 

The quoted uncertainties in equations ( |6.l[ ) and (3.2) 
represent the estimated overall errors in the model fits to 
the observed X-ray temperature function. The main source 
of uncertainty in our analysis is our modelling of the X-ray 
emitting intracluster medium as a homogeneous, isothermal 
gas in hydrostatic equilibrium (cf equation 2.2). There is 
some tentative evidence from ASCA data that the tem- 
perature in some clust ers may be declining at large radii 
(Markevitch et al. 1996). On the other hand, hydrodynamic 



For Qo = 1 this estimate is independent of the shape of 



simulations show that the isothermal assumption is a good 
approximation in the region where most of the X-rays are 
emitted, at least in the case where cooling fiows are ignored 
(Evrard 1990, Tsai, Katz & Bertschinger 1994, Navarro et 
al. 1995). Our modelling also required us to fix a value of 
the parameter /3, the ratio of specific galaxy kinetic energy 
to specific gas thermal energy, and we chose to set 13 = 1, 
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Figure 11. Predicted number counts of S-Z clusters with Y > 



10 "^h arcmin^ out to 



0.33 and 0.5 as a function of Qq. 



The dotted line corresponds to open models (Aq = 0) and the 
dashed line to flat models ( Aq = 1 — f!o)- The fractional difference 
between the models is smaller than in the corresponding plot for 
the temperature function because the counts are dominated by 
objects at low redshift where the models have been normalised to 
match the observed temperature function. 



consistent with the results of simulations. If bulk motions 
or magnetic stresses contribute to the support of the gas, 
a larger value of /3 would be appropriate. In this case, our 
results may be recast by scaling all temperatures inversely 
with (3 and this would lead to higher estimates of as- Finally, 
the estimated gas temperatures and our inferred values of erg 
may also be underestimated if small-scale inhomogeneities 
in the gas distribution affect the measured X-ray spectrum. 

The values of erg required to match the local abundance 
of clusters may be compared with measurements of the fluc- 
tuation ampUtude on larger scales, particularly with those 
inferred from the microwave background anisotropics in the 
COBE 2-year data. By assuming a specific shape for the 
fluctuation spectrum, the COBE results may be extrapo- 
lated to 8 h~^Mpc and, in principle, comparison with equa- 
tions ( |6.l[ ) and ( |3.2| ) provides a test of the assumed spectral 
shape. This test, however, cannot yet be made completely 
rigorous because the available anisotropy data do not dis- 
tinguish between contributions from scalar modes which de- 
termine the amplitude of the mass fluctuations, and con- 
tributions from tensor modes which produce gravitational 
waves. Furthermore, the asymptotic slope of the mass power 
spectrum, n, is poorly constrained by the anisotropy data 
and this introduces further uncertainty in the analysis of 
the COBE results (but not in the cluster abundance argu- 
ment). As an illustration, we compare, in Fig our esti- 
mates of (Tg with the COBE values for the simplest possible 



Figure 12. Sunyaev-Zel'dovich source counts as a function of Y 
for three different cosmologies. Solid lines correspond to f!o = 1, 
dotted lines to Qq = 0.3, and dashed lines to Qg = 0.3, Aq = 0.7. 
The corresponding 0.5mm fluxes in mjy are shown on the top 
axis. About 60% more clusters with Y > lO^^h arcmin^ are 
expected in the Oq = 0-3 models than in the Qq = 1 case. 



case in which n takes the standard value of unity and ten- 
sor contributions are neglected. The hashed regions in the 
Figure show COBE estimates for h = 0.8 and h — 0.5, 
taken from Liddle et al. (1995a) in the open case, and from 
Liddle et al. (1995b) in the flat case. The dashed and solid 
lines give our estimates of erg for the same two values of h, 
derived from the procedure described in Se ction 4, setting 
r — f2o/iexp(— f2b — Slb/f^o) ( Sugiyama 1995) with th e nucle- 
osynthesis value of f2b = 0.013/h (Copi et at 1995). Large 
values of h require low values of Qq. Thus, if ft = 0.8, Qo is 
required to be less than 0.2, in the flat case, and less than 
0.35 in the open case. If h > 0.5, values of Qo > 0.5 are ruled 
out by this test whether or not Aq = 0. For Qo = i, h — 0.27 
is required for consistency between COBE and the cluster 
abundances. Thus, as noted before by Efstathiou, Bond & 
White (1992) and White et al. (1993), the abundance of 
clusters in the standard CDM model {h = 0.5) is incompat- 
ible with the COBE fluctuations unless gravitational waves 
make a significant contribution to the measured microwave 
background anisotropics or there is a strong tilt in the pri- 
mordial spectrum (n < 0.8; Bond 1995). 

While the present day abundance of clusters may be 
used to determine the value of erg with only a weak depen- 
dence on the power spectrum of fluctuations, the evolution 
of the cluster abundance may be used to set constraints on 
f2o itself. These depend sensitively on the value of ag and, to 
a lesser extent, on the shape of the fluctuation spectrum. In 
this study, we have, for the most part, restricted attention to 
models with spectral shape parameter, F = 0.25, consistent 
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Figure 13. Comparison of erg estimated from the COBE 2-year 
data and from tlio abundance of rich galaxy clusters. The top 
panel corresponds to flat models and the bottom panel to open 
models. The upper hatched region shows the COBE estimates for 
h = 0.8 and the lower hatched region the estimates for h = 0.5 
as parameterised by Liddle et al. (1995a and b) with 12% er- 
rors in both cases. The dashed and solid lines show our esti- 
mates from the cluster abundance for h = 0.8 and h = 0.5 re- 
spectively, obtained through the procedure of Section 4, setting 
r = Qoh exp(— fib — Q^i/Qo), where Q^, = 0.013//i^ is the ratio 
of the mean baryon density to the critical density. The statistical 
errors in this case are less than 5%. 



with observations of the large-scale distribution of galaxies 
(Efstathiou et al. 1992). Since the Qo dependence of the evo- 
lutionary rate is so strong, the main effect of changing F is 
through its influence on the model normalization, ag. How- 
ever, as Figure |l^ shows, this is a weak effect. Thus, in prac- 
tice, this test of Oo is insensitive to small departures from 
our assumed value of F. Furthermore, for interesting values 
of Qo, the diagnostics we have considered depend only very 
weakly on Aq. Thus, statistical studies of clusters at interme- 
diate redshifts offer an excellent prospect for determining Qo 
without the complications arising from the uncertain shape 



of the galaxy power spectrum and the poorly understood 
relation between the distributions of galaxies and mass. 

Different observables may be used to characterize the 
evolution of the cluster population. Here we have consid- 
ered three properties for which observational data are likely 
to be obtained in the near future: the distribution of cluster 
masses. X-ray temperatures and S-Z decrements. The first 
may be derived from weak gravitational lensing studies; the 
second from existing and forthcoming X-ray surveys; and 
the last from ground-based observations and proposed space 
missions such as COBRA/SAMBA. Although the masses are 
perhaps the most difficult to measure, they provide a par- 
ticularly robust test since, apart from the Press-Schechter 
ansatz, the only model assumption is the form of the den- 
sity distribution of dark matter halos. In rough agreement 
with N-body simulations, we have taken this to be a singular 
isothermal sphere. In all cases, we have presented expected 
distributions at z = 0.33 and z = 0.5 and the predicted 
number counts of the largest clusters, both in space and in 
projection on the sky, as a function of redshift. 

In agreement with previous analyses, we find that even 
at redshifts as low as 0.33, the comoving abundance of high 
mass or high X-ray temperature clusters is radically different 
in universes with Qq ~ ^ and Qo = 0.3. For example, at this 
redshift we expect 20 times as many clusters with M > 3.5 x 
10^''h~^MQ and 5 times as many clusters with kT > 5 keV 
if Qo ~ 0.3 than if fio = 1- The corresponding factors at z — 
0.5 are ~ 65 and 15. When integrating these cluster number 
densities with respect to redshift the effects of the different 
volume elements are included. Thus, the expected numbers 
of clusters per square degree on the sky out to redshift 0.5 
diflfer by a factor of ~ 10 for both M > 3.5 x lO^h-^M© 
and kT > 5 keV in these two cosmologies. 

Whereas the mass and temperature functions only 
evolve dramatically if f2o = 1, the distribution of S -Z de cre- 
ments (as measured by the y-function of equation 2.10) de- 
clines rapidly with redshift in all cosmologies. This is mainly 
because, for the range of redshifts that we have considered, 
an object far away subtends a smaller angle on the sky than 
the same object placed nearby. Even so, at z — 0.33, the co- 
moving number density of clusters with Y — lO^'^h arcmin^ 
is 10 times higher if Qo = 0.3 than if Qo = 1. However, the 
overall abundances at this redshift are all very low because 
of the rapid decrease in the angular size of the clusters with 
redshift. The projected counts are dominated by low redshift 
clusters and, as a result, the expected excess of clusters with 
Y > 10~^h arcmin^ is only ~ 60% if Qo = 0.3. This excess 
becomes larger as the Y threshold is lowered and higher 
redshift clusters are included. 

To obtain an estimate of Qq using any of the tests pro- 
posed in this paper requires complete samples of intermedi- 
ate redshift clusters selected according to the statistic un- 
der consideration. Samples of S-Z clusters are likely to be 
y-limited in any case. For the other diagnostics, selecting 
by X-ray luminosity is probably an efficient method. Since 
X-ray luminosity correlates reasonably well with X-ray tem- 
peratu re (at least locally ) and, apparently also with lensing 
mass ( Smail et al. 1996 ), the sample completeness with re- 
spect to the relevant statistic may be determined a posteri- 
ori from these data. Incomplete datasets, on the other hand, 
may be used to set upper limits on Qo because the abundance 
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of clusters at intermediate redshift declines monotonically 
with increasing f2o- 

There are already several X-ray bright clus ters k nown at 
high redshift. For example, Luppino & Gioia ( 1995| ) report 
the dete ction of 6 such clu sters with z — 0.5 in the EMSS 
survey (Henry et al. 1992). Measurements of their X-ray 



temperatures would be extremely valuable. They might con- 
firm the indication from velocity dispersion measurements 
that these clusters are hot, suggesting that f2o is low, or 
they might show that these velocity dispersions are overes- 
timated as a result of contamination by projection effects. 
The lowest velocity dispersion reported for these clusters 
corresponds to a temperature of 6.6 keV. As an illustration, 
if we assume that the effective area surveyed to find these 
six clusters was 200 square degrees, then we would expect 
to find < 0.04 clusters above this temperature if f2o = 1; 
~ 6 ff r2o = 0.3 and Ao = 0.7; and ~ 30 ff Vlo = 0.3 and 
Ao = 0. However, if the coolest of these clusters actually has 
a temperature of only 3 keV, then we would have expected 
~ 30 examples even if f2o = 1- 

In summary, a statistical sample of X-ray clusters at 
intermediate redshifts with measured temperatures or S-Z 
decrements could place a strong constraint on the density 
parameter flo- 
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APPENDIX A: THE SPHERICAL COLLAPSE 
MODEL FOR + A = 1 



It can be shown (Peebles 1984) that the variation of the 
radius, r, of a uniform spherical overdensity with scale fac- 
tor a is described, for a flat cosmology with a cosmological 
constant, by 



dr \ ^ r ^ + ujr^ — k 



Vd^ 



(Al) 



•.da J a 1 + uja^ ' 

where a=(l + z)~^,Kisa constant which for overdensities 
takes positive values and 



(A2) 



Therefore a perturbed region will turnaround when 

_3 



-97 



Solving this cubic and requiring that a physically sensible 
root exists we find the smallest perturbation which will col- 
lapse has 



_ 3UJ3 
^min — 9 • 

2i 

The corresponding turnaround radius is 



(A4) 



(A5) 



Perturbations with larger values of k turnaround and col- 
lapse earlier. Solving the cubic (A3) we find the turnaround 
radius rta as a function of the density parameter ui and per- 
turbation amplitude k. The solution can be expressed as 



rta = — 2s3 cos 
where 

3 



/6» 27r\ 



s — 



and 9 satisfies 



0<.<^ 



(A6) 



(A7) 



(A8) 



In order to calculate the redshifts correspondin g to 
turnaround and collapse we separate the variables in (Al) 
and integrate. 



r 2 



/o {uir^ — Kr + 1) 
Defining 



-dr 



(tja^-fl): 



-dr, 



-do. 



(A9) 



(AlO) 



Iq (ujr-^ — Kr + 1) : 

which must be evaluated numerically, the expansion factor 
at turnaround, ata, and at collapse, a^, satisfy the implicit 
equations, 



and 



(wa3 -1-1)^ 



-da 



'o {uja^ + 1) 2 



-da. 



(All) 



(A12) 



where the collapse time is defined in the usual way to be 
that at which (Al) predicts collapse to a point singularity. 

these equations can be 



Using the substitution u — V^oa'' 
integrated analytically to give 



ata = 



and 



exTp{3^/ujl{uj , k)) — 1 



2v/ZJ 



exp{6\/uj 1(lo,k)) — 1 
2^ 



exp{-^I{iu,K)) (A13) 



exp{-2^I{u;,K.)) (A14) 



Hence we have Uc, or equivalently Zc, the redshift of collapse 
of a spherical overdensity, as a function of u and k. 

To calculate Ac, the density of virialised object in terms 
of critical dens ity, we make use of the expression from (La 



hav et al. 1991) which gives the ratio between the radius of 



the virialised sphere and its turnaround radius. Given the 
expression above for the turnaround radius and the expan- 
sion factor at collapse, Ac is given by 

Ac^n(ac) l — (A15) 



.(ac)(j)^ 



All that remains is to relate k to the quantity more 
often used to describe a perturbation, namely So, the linear 
theory overdensity extrapolated to the present. This in turn 
is related to the linear theory overdensity at collapse, 5c, by 



5c = 5oD{ac) 



(A16) 



where D{a) represents the linear theory growth factor nor- 
malised to unity at the present, D{ao) = 1. For the f2o+Ao — 
1, fio < 1 case 



D{a) = 



A{a{2ujy^'') 
A(ao(2L^)i/3) 



where 



A{x) 



(^^3-1-2) 



du 



(A17) 



(A18) 



(Peebles 1980). In the limit where a ^ 1 we find 

^000(2^)'/=^ 2. 
^^ 5A(ao(2a.)V3) " + Q(=^ )■ (^19) 

which can be compared to the corresponding expression 
3k 



-a + 0(a^ 



derived from (Al). Identifying terms we find 
ao(2u;)^^-' 



3A(ao(2Lj)i/3) 



5o. 



(A20) 



(A21) 



Now given the density parameter flo and the ampli tude of a 
density perturbation, 5o we can compute uj using (A2) and 
K using (A21). Then the collapse redshift can be obtained 
from (A6),(A10) and (A14 ). T he cor respo nding values of 5c 
and Ac are then given by (AW) and ( A15). The dashed lines 
in Figure ^ show 5c and Ac as functions of f2 at the time of 
collapse. 
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